arXiv: 1505.00817vl [astro-ph.SR] 4 May 2015 


Accepted for ApJ; 05/04/2015 


HST/COS Detection of the Spectrum of the Subdwarf 

Companion of KOI-81^ 

Rachel A. Matson^, Douglas R. Gies, Zhao Guo, Samuel N. Quinn 

Center for High Angular Resolution Astronomy and Department of Physics and Astronomy, 
Georgia State University, P. O. Box 5060, Atlanta, GA 30302-5060, USA; 
rmatson@chara.gsu.edu, gies@chara.gsu.edu, guo@chara.gsu.edu, guinn@astro.gsu.edu 

Lars A. Buchhave^, David W. Latham 

Harvard-Smithsonian Center for Astrophysics, Cambridge, MA 02138, USA; 
lbuchhave@cfa.harvard, edu, dlatham@cfa.harvard, edu 

Steve B. Howelh and Jason F. Rowe^ 

NASA Ames Research Center, P. O. Box 1, M/S 244-30, Moffett Field, CA 94035, USA; 

steve.b.howell@nasa.gov, Jason.Rowe@nasa.gov 

ABSTRACT 

KOI-81 is a totally eclipsing binary discovered by the Kepler mission that 
consists of a rapidly rotating B-type star and a small, hot companion. The 
system was forged through large scale mass transfer that stripped the mass donor 
of its envelope and spun up the mass gainer star. We present an analysis of UV 
spectra of KOI-81 that were obtained with the Gosmic Origins Spectrograph on 
the Hubble Space Telescope that reveal for the first time the spectral features of 
the faint, hot companion. We present a double-lined spectroscopic orbit for the 
system that yields mass estimates of 2.92Mq and 0.19Mq for the B-star and hot 
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subdwarf, respectively. We used a Doppler tomography algorithm to reconstruct 
the UV spectra of the components, and a comparison of the reconstructed and 
model spectra yields effective temperatures of 12 kK and 19 - 27 kK for the B- 
star and hot companion, respectively. The B-star is pulsating, and we identihed 
a number of peaks in the Fourier transform of the light curve, including one 
that may indicate an equatorial rotation period of 11.5 hours. The B-star has 
an equatorial velocity that is 74% of the critical velocity where centrifugal and 
gravitational accelerations balance at the equator, and we £t the transit light 
curve by calculating a rotationally distorted model for the photosphere of the 
B-star. 

Subject headings: stars: individual (KOI-81) — stars: binaries: spectroscopic — 
stars: evolution — stars: subdwarfs 


1. Introduction 

The NASA Kepler mission of precise and long time span photometry has led to many 
surprising discoveries. Among the initial set of remarkable targets, Rowe et ah (2010) re¬ 
ported on observations of two transiting systems dubbed KOI-74 (P = 5.2 d) and KOI-81 
(P = 23.9 d) that display light curves with minima that were deeper during occultation 
than during transit, implying that the planetary size companions are hotter than their A- or 
B-type host stars. Rowe et ah (2010) proposed that these companions were actually very low 
mass white dwarf (WD) stars, the remnants of more massive progenitors in an interacting 
binary. The light curves of both systems were further analyzed by van Kerkwijk et al. (2010) 
who found that both appear more luminous during orbital phases of the bright star approach 
due to relativistic Doppler boosting (“beaming binaries”; Zucker et al. 2007). Van Kerkwijk 
et al. used the radial velocity semiamplitude implied by Doppler boosting to estimate that 
the white dwarf stars had masses of ~ 0.2 — O.SMq. 

The Aep/er observations eventually led to the discovery of other similar eclipsing binary 
systems. Carter et al. (2011) found that KIC 10657664 consists of a hot WD orbiting an 
A-type star (P = 3.3 d), and Breton et al. (2012) discovered a WD and F-star binary system, 
KOI-1224 (P = 2.7 d). Two more systems orbiting A-type stars were recently discovered by 
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Rappaport et al. (2015), KIC 9164561 (P = 1.3 d) and KIC 10727668 (P = 2.3 d). Maxted et 
al. (2011, 2014) reported on the detection of 18 similar short-period systems (P = 0.7 —2.2 d) 
in photometric time series from the ground-based Wide Angle Search for Planets survey 
program. Photometric detection of longer period systems is less probable because the orbital 
inclination must be extremely close to 90° for transits to occur. However, longer period 
systems can be found through extensive, time series, radial velocity measurements from 
spectroscopy. For example, Gies et al. (2008) used radial velocity measurements to show 
that the nearby star Regulus is a spectroscopic binary (P = 40.1 d) consisting of a B-star 
with a probable WD companion. These hot companions of main sequence stars are probably 
related to the subdwarf companions of rapidly rotating Be stars that are detected through 
ultraviolet spectroscopy (Peters et al. 2013). 

Close binaries are common among intermediate mass stars, and many of these will 
experience large scale mass transfer (Rappaport et al. 2009; van Kerkwijk et al. 2010; Di 
Stefano 2011; Clausen et al. 2012). As the initially more massive star grows in radius, this 
donor star will begin Roche lobe overflow (RLOF) and transfer mass and angular momentum 
to the companion (the mass gainer). The orbit will shrink until the masses are equal, and 
if contact can be avoided, then additional mass transfer will cause an expansion of the 
orbit that ceases once the donor has lost its envelope. The resulting system will consist of 
a stripped-down and hot donor star in orbit around a rapidly rotating and more massive 
gainer star. 

Detecting binaries in this stage of evolution is difficult because the small donor stars 
are relatively faint and lost in the glare of the brighter companions. Furthermore, the 
donors are now low mass objects that create only small reflex orbital motions in the gainer 
stars, so detection through spectroscopy is challenging. Thus, the Kepler discoveries offer 
us a remarkable opportunity to investigate the properties of stars that are known examples 
of this post-mass transfer state. The hot companions will contribute a larger fraction of 
the total flux at shorter wavelengths, so a direct search for the flux and spectral features 
associated with the hot companion is best done in the ultraviolet. Indeed, it was through 
UV spectroscopy from the International Ultraviolet Explorer satellite (Thaller et al. 1995) 
and from the Hubble Space Telescope (Gies et al. 1998) that the hot subdwarf companion of 
the Be star 0 Per was first detected. 

Here we report on an UV spectroscopic investigation of KOI-81 (KIC 8823868; TYC 3556- 
3094-1; 2MASS J19350857-I-4501065) made possible with the HST Cosmic Origins Spectro¬ 
graph (COS) that has revealed the spectral features of the hot companion for the first time. 
We describe the COS observations and supporting ground-based spectroscopy in §2. We 
present radial velocity measurements and a double-lined orbital solution in §3 to obtain 
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mass estimates. We then use the derived radial velocity curves to perform a Doppler tomo¬ 
graphic reconstruction of the component spectra, and we compare the reconstructed spectra 
to model spectra to derive effective temperatures and projected rotational velocities (§4). 
The Kepler light curve outside eclipses shows evidence of pulsational and rotational fre¬ 
quency signals that we discuss in §5. The transit light curve is analyzed in §6 using a model 
for the rotationally distorted B-star. We summarize and consider some consequences of the 
results in 57. 


2. Spectroscopic Observations 

The Cosmic Origins Spectrograph (COS) is a high dispersion instrument designed to 
record the UV spectra of point sources (Osterman et al. 2011; Green et ah 2012). The 
HST/COS observations of KOI-81 were obtained over hve visits between 2011 June and 
2011 October. These were scheduled so that two occurred during each of the quadrature 
phases near the Doppler shift extrema, and the hfth observation was made during the hot 
star occultation phase in order to isolate the flux of the B-star alone. The UV spectra were 
made with the G130M grating to record the spectrum over the range from 1150 to 1450 A 
with a spectral resolving power of R = A/AA = 18000. There are two COS detectors that are 
separated by a small gap, so the spectra were made at slightly different central wavelengths 
in order to hll in the missing flux: 1300 A (four exposures of 447 s), 1309 A (three exposures 
of 399 s), and 1318 A (three exposures of 399 s). This sequence required an allocation of two 
orbits for each visit. The observations were processed with the standard COS pipeline to 
create wavelength and flux calibrated spectra as xld.fits hies for each central wavelength 
arrangement (Massa et al. 2013). These ten sub-exposures were subsequently merged onto 
a single barycentric wavelength grid using the IDL procedure coadd_xld.prq^ (Danforth et 
al. 2010). We created a list of the sharp interstellar lines in the spectrum, and we cross- 
correlated each of these spectral regions with those in the average spectrum in order to make 
small corrections to the wavelength calibration. Then the interstellar lines were removed 
in each spectrum by linear interpolation across their prohles. Finally, all hve spectra were 
transformed to a uniform grid with a log A pixel spacing equivalent to a Doppler shift step 
size of 2.60 km s“^ over the range from 1150 to 1440 A. The coadded spectra have a signal- 
to-noise ratio of S/N = 110 in the central, best exposed parts. 

We also obtained three sets of complementary ground-based spectra of KOI-81. The 
hrst set consists of 19 high dispersion spectra made with Tillinghast Rehector Echelle Spec- 
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trograph (TRE^ mounted on the 1.5 m Tillinghast telescope at the Fred Lawrence Whipple 
Observatory at Mount Hopkins, Arizona. These spectra cover the optical range from 3850 
to 9100 A with a resolving power oi R = A/AA = 48000. The spectra were processed and 
rectihed to intensity versus wavelength using standard procedures (Buchhave et ah 2010), 
and th^ are available at the Kepler Community Follow-up Observing Program (CFOP) 
websitq^- A second set of six moderate resolution spectra were made in 2010 and 2012 with 
the Kitt Peak National Observatory (KPNO) 4 m Mayall telescope and the Ritchey-Chretien 
(RC) Focus Spectrograph. These were made with the BL 380 grating (1200 grooves mm“^) 
to record the spectrum between 3950 to 4600 A with a resolving power of i? = 6300. Finally 
a third set of hve, lower resolution, and flux calibrated spectra were obtained in 2010 with 
the Mayall telescope and RC spectrograph using the KPC-22B grating (632 grooves mm“^) 
to cover the region from 3577 to 5058 A with a resolving power of i? = 2500. This third set 
is also available at the CFOP website. 


3. Radial Velocities and Orbital Elements 

Our first goal was to measure the orbital motion of the primary star using the COS UV 
spectra. Direct inspection of the spectra showed that the main features were very broadened 
and blended due to large rotational broadening (§4). Consequently, we decided to measure 
the individual spectrum velocities by cross-correlating them with a spectral template. We 
formed cross-correlation functions (CCFs) using the hot star occupation phase spectrum 
(made on HJD 2455775.9626) as the CCF template. The calculation was made using only 
the regions from 1270 to 1300 A and from 1314 to 1437 A (see Fig. 4 below) in order to 
avoid the low wavelength region where the secondary’s flux becomes larger and to remove 
the Si II AA1260,1264,1304,1309 features that may be affected by incomplete removal of 
the interstellar components. The resulting CCFs are extremely broad, and we found the 
center of each CCF by measuring the wing bisector position by convolving the CCF wings 
with oppositely signed Gaussian functions (Shatter et al. 1986). We similarly measured the 
center of the CCF of the hot star occupation phase spectrum with a model spectrum for the 
B-star (§4) from the UVBLUE library (Rodriguez-Merino et al. 2005), and this offset was 
added to the relative velocities to transform them to an absolute scale. The resulting radial 
velocities are collected in Table 1 that lists a leading P or S for primary B-star or secondary 
hot star (see below), the heliocentric Julian date of mid-observation, the orbital phase based 
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upon the ifep/er solution for the time of mid-transit (B-star superior conjunction), the radial 
velocity, the difference between the radial velocity and the systemic velocity for the specihc 
observation set (see below), the uncertainty a, the observed minus calculated (0—C) velocity 
residual from the £t (§4), and the source (HST/COS in this case). Both pairs of qnadrature 
observations were separated in time by only a few hours and are expected to have almost 
the same orbital velocity, so we used the absolute differences of the measured velocities of 
these pairs to estimate \/2a. 

We realized at the outset that finding and measnring the spectral lines of the companion 
would be challenging because of its relative faintness. However, detection is favored at 
the shorter wavelength part of the UV spectrum because the hotter companion contribntes 
relatively more flux there. An inspection of the spectra in the region with A < 1200 A did 
indeed hint at the presence of a velocity variable, narrow-lined component. We needed to 
isolate the flux contribution from the bright B-star in this region in order to remove its flux 
from each spectrum and reveal the secondary’s spectral lines in the flux difference. This was 
accomplished using the hot star occultation phase spectrum to represent the flux of the B- 
star alone. This spectrnm was smoothed by convolution with a Gaussian function of FHWM 
= 133 km s“^ in order to increase the S/N ratio without unduly altering the spectral shape. 
Then the smoothed version was shifted in velocity according to the predicted orbital motion 
of the primary star, and a difference spectrum was formed for each of the four quadrature 
phase observations. Renormalized and averaged versions of these qnadratnre phase spectra 
are illustrated in Figure 1, and these show a host of weak and narrow lines that are Doppler 
shifted as expected for the orbital motion of the hot secondary. The flux scaling for the hot 
companion is not accurate because of the rudimentary means of the B-star’s flux removal, 
but it is sufficient to reveal the lines of the hot component so that its radial velocity can 
be measured. We did so by calculating the CCF of each qnadrature spectrum using a hot 
model spectrum (§4) for the template, and the center of the narrow signal in the resulting 
CCF was estimated by fitting a parabola to its peak (with uncertainties determined using 
the method of Zncker 2003). The derived radial velocities of the secondary are listed at 
bottom of Table 1 in rows with a leading colnmn marked S. No measnrement is given for the 
occultation phase spectrum because the flux of the hot component is totally blocked then 
(Rowe et ah 2010). 

We also measnred radial velocities for the primary B-star using the higher resolntion, 
ground-based spectra. The 19 TRES spectra have a S/N that is too low for measurement 
of the weak and broadened He and metallic lines, but they are sufficient to measure the 
Doppler shifts in the strong and broad, H Balmer lines. We formed CCFs of each of Ha, 



7 


H/?, H 7 , H(5, and He with model template spectra from the BLUEREDIj grid (Bertone et al. 
2008). The radial velocity was estimated by htting a parabola to the peak of the CCF for 
each line. The qnality of the measnrement varies depending on the position of the featnre 
relative to the echelle blaze fnnction (best for Ha and H 7 which fall close to the blaze 
maximnm), so we formed weights for each of the Balmer lines proportional to the inverse 
sqnare of the standard deviation of the velocity measnrements. We then determined the 
weighted mean and standard deviation of the mean from the set of five line measnrements 
for each observation, and these are listed in Table 1 (noted by TRES in the final colnmn). 

We followed a similar procednre to measnre radial velocities for the six KPNO moderate 
resolntion spectra. We calcnlated CCFs nsing a model B-star template (from the UVBLUE 
grid) and inclnding all the lines in the region from 4050 to 4520 A (dominated by H 7 and 
H5), and we fit a parabola to the CCF peak to estimate radial velocities. These are identihed 
in Table 1 by KPNO in the final colnmn. 

We fit orbital elements for each of the primary and secondary velocity sets nsing the 
program described by Morbey & Brosterhus (1974). We fixed the orbital period and epoch 
of mid-transit to those derived by the Kepler project (dated 2014 November 24) that are 
posted at the CFOP website. Because the transit and occultation light curves have the 
same duration and are separated by half of the orbital period, we assumed that the orbit 
is circular. We set the fitting weights to the inverse square of the uncertainties. Because 
the COS, TRES, and KPNO measurements are derived from differing spectral ranges and 
instruments, we might expect that there will be systematic differences in the velocities for 
each set. The upper panel of Figure 2 shows the measured radial velocities from Table 1, 
using different symbols for each set of spectra, as well as preliminary fits of each set, and 
there do indeed appear to be constant offsets between the sets. We found independent 
estimates for the systemic velocity 7 * of each set by iterating between global (all data) and 
individual sets of velocities for the orbital fits. This was done by fixing the semiamplitude 
Ki from the global £t to find estimates of 7 * for solutions to each set. Then a new global 
fit was made using velocity differences (subtracting the 7 ^ value associated with each set) to 
hnd a revised Ki. This procedure converged after a few iterations to the final adopted fit 
presented in Table 2. The 7 j-corrected velocities (listed under column W ~ 7 i in Table 1) 
and 2a range of the fitted velocity of the primary are shown in the lower panel of Figure 2, 
and the combined velocity curves of the primary and secondary are illustrated in Figure 3. 
We also made a fit of the 7 i-corrected velocities of the primary using a Markov Chain Monte 
Carlo algorithm that led to almost the same result, Ki = 6 . 75 ^q 44 km s“^. However, we 
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conservatively adopt the somewhat larger uncertainty estimates from the non-linear least 
squares program in what follows. Van Kerkwijk et ah (2010) estimated ~ 7 km s“^ from 
the Doppler boosting in the light curve of KOI-81, and their result is verihed through our 
direct Doppler shift measurement of Ki = 6.7 ± 0.7 km s“^. 

The orbital inclination is very close io i = 90°. We derive a value from the transit light 
curve (§6) oi i = 88?97 ± 0?04, which is intermediate between the estimates from Rowe et ah 
(2010) oii = 88?2±0?3 and from the i7ep/er project posted at CFOP oii = 89?95. Thus, we 
can use our estimate to derive the physical masses from the Mi sin^ i, M 2 sin^ i products, and 
these are given in Table 3. Furthermore, the average density p of the B-star can be directly 
estimated from the transit light curve (provided the star is spherical; see §6), and the Kepler 
project hnds p = 0.280 ± 0.005 g cm“^ for KOI-81 (reported at the CFOP website). We 
used this value together with the mass estimate Mi to arrive at the radius Ri reported in 
Table 3. Finally, the CFOP website gives the ratio of the radii derived from the transit light 
curve, R 2 /Ri = 0.03725 ± 0.00026, and we used this ratio to find R 2 (Table 3). Table 3 also 
presents the gravitational acceleration logg derived from the mass and radius information. 
We must examine the spectrum of the system to derive the remaining stellar parameters. 


4. Tomographic Reconstruction of the UV Spectra 

We used a Doppler tomography algorithm (Bagnuolo et al. 1994) to extract the indi¬ 
vidual UV spectra of the primary and secondary stars. This is an iterative scheme that uses 
estimates of the orbital Doppler shifts of each component and their flux ratio to derive re¬ 
constructed spectra for both stars. We initially assumed featureless continua as the starting 
approximation for the spectra of both stars, but after comparison with model spectra, we ran 
the algorithm again using the models as starting values, and this choice helped to limit the 
continuum wander in the resulting spectral reconstructions. We present in Figures 4 and 5 
the reconstructed UV spectra derived from the four COS spectra obtained at the quadrature 
phases. Figure 4 also shows the excellent agreement between the reconstructed spectrum 
of the primary and the hot star occultation phase spectrum that represents the flux of the 
primary alone. Figures 4 and 5 indicate some of the principal lines as well as the locations 
of the interstellar lines, where their removal by interpolation may have interfered with the 
accurate reconstruction of the spectra. All the spectra in Figures 4 and 5 are smoothed by 
convolution with a boxcar function of width 133 km s“^ in order to reduce the noise and 
facilitate intercomparison of the line features. 

We compared the reconstructed spectra with model spectra from the UVBLUE grid of 
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high resolution spectralj calculated by Rodriguez-Merino et ah (2005) that are based upon 
the ATLAS9 model atmosphere code and SYNTHE radiative transfer code developed by 
R. L. Kurucz. We used their solar metallicity models that incorporate a microturbulent 
velocity of 2 km s“^. We made simple bilinear interpolations of flux in the (Tefj, loggf) plane 
to derive the model spectra. We adopted log 5^1 = 4.13 for the primary (Table 3), but set 
log 5^2 = 5.0 for the secondary because this value is the largest available in the UVBLUE 
grid. The model spectra were rebinned onto the observed wavelength grid and then convolved 
with the instrumental broadening function (from the COS line spread function!^ for a central 
wavelength of 1300 A) and with a rotational broadening function (using linear limb darkening 
coefficients from Wade Sz Rucinski 1985). 

We first compared the reconstructed and model spectra to estimate the projected rota¬ 
tional velocity V sini of each star. This was done by forming a goodness-of-£t statistic 
between the observed and model spectra over a test grid of V sin i values, and then finding 
the V sini corresponding to the minimum This was repeated for a series of wavelength 
regions that contained well defined absorption lines or line blends, and the mean and stan¬ 
dard deviation of the derived V sini for the primary is presented in Table 3. The primary 
B-star is indeed a very rapidly rotating star, a property noted first by van Kerkwijk et ah 
(2010). The lines of the hot secondary companion, on the other hand, appear very sharp and 
any rotational broadening is unresolved in the COS spectra, so we present only an upper 
limit in Table 3. 


Next we needed to adjust the model spectra for the variable interstellar extinction 
across the COS wavelength band. We derived the interstellar reddening by comparing the 
available observed fluxes with a model flux distribution for the binary transformed using 
the extinction law presented by Fitzpatrick (1999). The observed fluxes are shown in a 
spectral energy distribution plot in Figure 6 . These include rebinned COS measurements, 
a GALEN0 NUV measurement!^, rebinned KPNO low dispersion spectrophotometry, and 
2MASS and WISE photometric fluxes. The model fluxes were a sum of Kurucz models 
for the primary and secondary that were attenuated according to the extinction law for a 
ratio of total-to-selective extinction of R = 3.1. The £t shown in Figure 6 was made with 
a reddening of E{B — V) = 0.169 ± 0.008 mag and a limb darkened angular diameter of 


®http://www.inaoep.mx/'^modelos/uvblue/go.htinl 
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^^KOI-81 appears as two sources in GALEX, but we simply summed the fluxes of what appears to be the 
trailed image of a single object. 



- 10 - 


the primary of 6^ = 0.0182 ± 0.0005 milliarcsec. The former agrees with prior estimates 
{E{B — V) = 0.175 and 0.193 mag in the Kepler Input Catalog and the GALEX Catalog, 
respectively), and the latter indicates a distance of 1.25 kpc (using Ri/Rq from Table 3). 
The extinction in the COS band was calculated for the derived value of E{B — V), and 
each model was multiplied by this extinction curve to account for the greater interstellar 
attenuation at lower wavelength. 

With the gravity, projected rotational velocity, and interstellar extinction set, we then 
calculated model spectra over a grid of effective temperature and determined a goodness- 
of-£t as a function of T^s- The best fit temperatures and their estimated uncertainties from 
this comparison of the reconstructed and model spectra are listed in Table 3. These un¬ 
certainties do not account for possible differences from the assumed microturbulent velocity 
and abundance. Nevertheless, these models offer a framework to interpret the appearance 
of the UV spectra. 

Figure 4 shows the derived model spectrum offset above the reconstructed spectrum of 
the primary star. The overall agreement is very good in both the appearance of the lines 
and line blends. The C II AA1335,1336 feature is weaker than predicted, but this may result 
from our means of removal of interstellar line components in this feature. The flux in the 
optical wavelength range is totally dominated by the light from the B-star primary (§6), so 
we did not attempt Doppler tomography reconstructions of the optical spectrum. However, 
we show in Figure 7 two optical lines of special interest. The most temperature sensitive 
feature in the optical spectrum at Tefj ~ 10 kK is the Ca II A3933 K line that grows rapidly in 
strength with decreasing temperature (Gray & Corbally 2009). The average TRES spectrum 
of this feature is shown in Figure 7 {left) together with model spectra for Rs = 11.7 kK and 
10.6 kK, and a comparison with the observed K-line suggests that the effective temperature 
may be somewhat lower than that derived from the ultraviolet spectrum (but still within 
the uncertainties). On the other hand, the difference may result from gravity darkening in 
the rapidly rotating B-star that creates a cooler equatorial region (§6) that would promote 
the appearance of lines favored in cooler atmospheres. Figure 7 also illustrates the average 
observed and model profiles of the gravity sensitive Hy line {right), and the observed profile 
has a shape that is approximately consistent with predictions for the adopted gravity of the 
primary. 

The reconstructed secondary spectrum shown in Figure 5 is quite noisy because this 
component is relatively faint (see below), but there is satisfactory agreement between the 
reconstructed and model spectra in the vicinity of some features. For example, the recon¬ 
structed and model spectra of the secondary appear similar in the lines of C III A1176, 
Si II AA1193,1194, Si III AA1295,1297,1299, and Si IV AA1394,1403. The latter are only 
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seen in the spectrum of the hotter secondary. However, the agreement is less satisfactory 
in the regions where interstellar lines were removed (near 1229, 1251, 1303, 1335, and 1347 
A). The model spectra were formed using UVBLUE models with loggf = 5.0, the largest 
value in the grid, which is smaller than the estimated gravity, logg = 5.8. This has two 
important consequences. First, the model Lya line is narrower than observed because the 
linear Stark broadening associated with the lower gravity model is insufficient to match 
the observed broadening. Second, the effective temperature estimation is based upon the 
relative strengths of transitions corresponding to different ionization states, and according 
to the Saha equation, line formation in a denser medium (at higher gravity) will shift the 
ionization balance to less ionized states. Consequently, a good £t of the spectrum would 
also be possible with a higher gravity and higher temperature model, so our derived effective 
temperature should be regarded as a lower limit. We also show in Figure 5 an example 
spectrum of the B-subdwarf CPD—64°481 that was made with the Space Telescope Imaging 
Spectrograph and E140H grating. O’Toole & Heber (2006) used this and other spectra to 
estimate Tes = 27.5 kK and logg = 5.6, a value of gravity that is closer to that of the sec¬ 
ondary in KOI-81. Comparing these similar gravity stellar spectra, we see that transitions 
of Si II are almost absent in the spectrum of CPD—64°481, whereas they are quite strong 
in the spectrum of the KOI-81 secondary, so the KOI-81 secondary must be cooler than 
CPD—64°481. Thus, we can place lower and upper limits on the effective temperature of 
the secondary of 19 < < 27 kK. 

All the spectra in Figures 4 and 5 were normalized in flux by dividing by the mean 
flux over the range between 1340 and 1400 A. We set the monochromatic flux ratio from 
this normalization to T2(1370)/Fi(1370) = 0.05 ± 0.02 based upon a comparison of the line 
depths in the reconstructed and model spectra of the secondary. This comparison was made 
using relative line depths in the vicinity of the lines that matched reasonably well, because 
there is too much wander in the continuum of the reconstructed spectrum of the secondary 
for a reliable global fit. We caution that the model line depths are sensitive to the adopted 
values of temperature, gravity, and especially microturbulence, so the flux ratio may need 
revision if different assumptions are made. 


5. Non-orbital Frequencies in the Light Curve 

Both Rowe et ah (2010) and van Kerkwijk et ah (2010) noted that the light curve of KOI- 
81 shows evidence of pulsations, and it is interesting to consider the pulsational frequencies 
and their relation to the rotational frequency of the B-star primary. We calculated the Fourier 
transform amplitude of the detrended, long cadence Kepler observations from quarters 0 
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to 17 using the package Period04 (Lenz & Breger 2005). The portions of the light curve 
covering the transit and occnltation of the secondary were removed from the time series in 
order to focus on non-orbital timescales of variation. We adopted an empirical noise level by 
smoothing the envelope of the residnal spectrnm after prewhitening all signihcant frequencies 
with S/N > 3 (except for the broad rotational featnre; see below). The nncertainties were 
calculated following Kallinger et al. (2008). The Fonrier spectrnm is dominated by a very 
strong signal with a freqnency of /i = 0.722974 cycles d“^ (period of 1.38318 d), so we 
removed this signal by pre-whitening to nncover the remaining periodic signals that are 
plotted in Figure 8. There are a number of strong signals that appear well above the noise 
level, and the freqnencies, amplitndes, sinusoidal phases (relative to the epoch of central 
transit), and peak signal-to-noise ratio are listed in Table 4. The hnal colnmn of Table 4 
lists several nnmerical relations among these frequencies, and identihes one freqnency fu 
that corresponds to the ellipsoidal (tidal) variation with half the orbital period (discussed 
by van Kerkwijk et al. 2010). 

The dominant periodicity (1.38 d) is nnrelated to the orbital period, and it probably 
represents a strong and long-lived pnlsational mode. This and the other low freqnency 
pulsation signals may be of two possible kinds. First, the B-star primary of KOI-81 has 
a temperatnre and radius that are similar to those of the slowly pulsating B-type (SPB) 
stars (Pamyatnykh 1999), and these stars display relatively long period gf-mode oscillations. 
Second, the B-star is a rapid rotator, and Townsend (2005) and Savonije (2013) argue that 
rapidly rotating, late-type B-stars can experience retrograde mixed modes of low azimnthal 
order m. The rotational frequency of the B-star is probably 2 cycles d“^ (see below), so the 
smaller freqnency of the dominant signal is consistent with retrograde nonradial pnlsation. 

There is a broad distribntion of peaks just above 2 cycles d“^, and we show an enlarged 
version of the Fonrier amplitude in this vicinity in the top panel of Fignre 9. A wide distri¬ 
bntion appears aronnd 2.04 cycles d“^ that is accompanied by a strong peak at /s = 2.08287 
cycles d“^. The lower panel shows the residnal peaks after pre-whitening and removal of 
the strong /s signal, and this reveals the presence of several other signihcant peaks near /s. 
This kind of broad feature with a stronger single peak at slightly higher freqnency has been 
detected by Balona (2013, 2014) in the Kepler light curves of some 19% of A-type stars. 
Balona argues persuasively that this feature probably corresponds to the stellar rotational 
freqnency. In his interpretation, the single strong peak corresponds to the eqnatorial rota¬ 
tional freqnency and the wider peak samples the rotational freqnencies at different latitndes 
in stars with differential rotation. Thus, following this line of argnment, we may tentatively 
identify /s as the equatorial rotational freqnency of the B-star in KOI-81, and thus, the ro¬ 
tational period is 0.48 d at the eqnator. Balona (2014) considers several explanations for the 
origin of the photometric variation inclnding pulsation, rotational modulation by starspots. 
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and tidal variations indnced by a co-orbiting exoplanet. We suspect that in the case of 
KOI-81, any co-orbiting planet would have a short period and an orbital plane similar to 
that of the stellar companion, so that we might expect to observe a transit signal in the /s 
folded light curve, but instead the folded light curve is approximately sinusoidal in shape. 
We speculate that the rotational signals in the light curve of KOI-81 and similar stars may 
result from long-lived vortices (Kitchatinov & Rudiger 2009) that develop in the outer atmo¬ 
spheres of rotating stars due to differential rotation (similarly to the spots in the atmosphere 
of Jupiter). 


6. Transit Light Curve 

The B-star we observe was spun up during the mass transfer process to produce a very 
rapidly rotating star with a rotationally broadened spectrum. It is important to consider 
how this rapid rotation influences our interpretation of the light curve. We expect that the 
spun up star will have a rotational axis parallel to the orbital angular momentum vector, 
so that the star’s spin axis also has an inclination of i ~ 90°. We argued in §5 that the 
photometric signal /s is the rotational frequency of the B-star. Then we may estimate the 
star’s equatorial radius from Ri equator = (Ksini)/(27r/5sini) = (2.81 ± O.O5)i?0. This is 
somewhat larger than what we derived from the mass and mean stellar density (Table 3), 
but is not unexpected for a rotationally distorted star in which the equatorial radius will be 
larger (and the polar radius smaller) than the mean radius. 

We can model the predicted appearance of the B-star based upon our estimates of stellar 
mass, equatorial radius, rotational period, and average effective temperature. We created 
a model image of the specihc intensity at a wavelength of 6430 A, which is the centroid 
of the Kepler instrument response function, using the same methods applied in a study of 
the rapidly rotating B-star Regains (McAlister et al. 2005). The star is assumed to have 
a shape that follows the Roche potential for rotation about a point mass, and each surface 
element has a specihc intensity dehned by limb and gravity darkening. The specihc intensity 
of a surface element is set by interpolation in a set of Kurucz model values dehned by the 
local temperature, gravity, and orientation of the surface normal to the line of sight. We 
found values of the polar radius and temperature that led to our derived equatorial radius 
and surface average temperature. An image of the rotationally distorted star appears in 
Figure 10. The resulting model has a ratio of equatorial to critical velocity (where the 
gravitational and centrifugal accelerations are equal) of 0.74, and the radius varies from 
2.3Oi?0 at the poles to 2.81i?0 at the equator. The temperature in this model varies between 
13.8 kK and 10.8 kK from pole to equator, and gravity likewise varies from logg = 4.18 to 
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3.75. The flux weighted, disk integrated value of gravity is \ogg = 3.9, and this appears to 
be approximately consistent with the appearance of the Hy profile (Fig. 7). 

We used this model image to calculate transit light curves. We assumed that the 
apparent transit path follows a trajectory parallel to the stellar equator with a relative 
velocity given by 27ra/P. During phases when the disk of the companion is seen projected 
against the background B-star, we took the flux removed as the product of the companion 
area and the specific intensity at the projected position of the center of companion (Mandel 
& Agol 2002; Barnes 2009). For the phases between first and second contact, we assumed 
that the companion occults a locally linear stellar limb of the B-star with a slope given 
by the Roche model shown in Figure 10, and the specific intensity was estimated at test 
positions perpendicular to the local limb to make a numerical calculation of the occulted 
flux. This constrained model has only two free parameters: the ratio of secondary to primary 
polar radius R 2 /R 1 polar and the ratio zq/Ri polar of the smallest distance zo from the transit 
trajectory to the center of the B-star relative to the primary polar radius. The value of 
R 2 /R 1 polar sets the amount of light removed and hence the depth of the transit light curve, 
while the value of Zq/Ri polar sets the duration of the transit (longer for transits near the 
equator with smaller zq). 

We show in Figure 11 the occultation and transit light curves of KOI-81 from nine 
months of Kepler short cadence observations. These plots show the photometric fluxes 
normalized to unity outside of eclipses that we calculated by rebinning and averaging all the 
available measurements in orbital phase bins equivalent to 3 minutes duration. The Kepler 
photometry was detrended and pre-whitened for the primary oscillation frequency /i (§5) 
before rebinning. The upper plot of the total occultation of the hot companion by the B-star 
indicates that the companion contributes a flux fraction of F 2 /F 1 = 0.00501 ± 0.00006 in 
the Kepler band-pass, and we renormalized the model for this extra flux in calculating the 
model transit curves. The model curves were also convolved with a temporal box function to 
represent the bin size applied to the observations. Our best fit (shown in Fig. 10 and Fig. 11) 
was obtained with R 2 /R 1 polar = 0.0426 ± 0.0003 and Zo/Ri polar = 0.396 ± 0.011 (equivalent 
to an orbital inclination of i = 88?97 ± 0?04). The individual transits display asymmetries 
that appear to be related to the pulsational phases at the times of transit, and because the 
mean transit curve shown in Figure 11 represents data from only 12 transits, we do not 
attribute any significance to the trends in the residuals from the fitted transit light curve. 
We caution that the quoted uncertainties in the fitting parameters do not account for the 
range in the possible values for the rotation rate we adopted. Nevertheless, it is interesting 
to note that the derived radius of the companion R 2 /Rq = 0.0979 ± 0.0019 is close to what 
we list in Table 3 based upon the spherical approximation, mean density, and radius ratio 
from the Kepler project, and the primary radius given in Table 3 falls comfortably between 
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the polar and equatorial radii of the rotating model. 


In §4 we compared the observed and model spectral line depths of the secondary to 
estimate that the total flux ratio at 1370 A is F 2 /F 1 = 0.05±0.02, and using the temperatures 
we derived, we predict that this flux ratio will vary from F 2 /F 1 = 0.14 at 1150 A to F 2 /F 1 = 
0.0076 in the i7ep/eroptical band. The observed flux ratio is approximately related to radius 
ratio by 


Fi h \Ri) 


( 1 ) 


where / 2//1 is the ratio of monochromatic flux per unit area. According to the Kurucz flux 
models for our stellar parameters, / 2//1 = 14.5 ± 0.5 at 1370 A, so the observed flux ratio 
would then imply a radius ratio of R 2 /R 1 = 0.059 ± 0.012. Both the flux and radius ratio 
estimates from spectroscopy are somewhat larger than the corresponding estimates from the 
Kepler transit analysis, but the differences are not surprising given the large uncertainties 
associated with the flux ratio from spectroscopy. 


We can also use equation 1 to infer the temperature ratio from the observed flux ratio 
in the Kepler band. The Kurucz flux models vary as / oc with a = 1.58 for logg = 5, 
A = 6430 A, and Teg in the range 10 to 25 kK. Thus, if we set the radius of the primary as 

{Rl equator + Rl polar)/2 = l.lli?i polar, then We CaU USe the values of F 2 /F 1 and R 2 /R 1 polar 
from above to find a temperature ratio of T2 es/Ti eg = 2.2. This is broadly consistent with 
the temperature estimates from spectroscopy (§4), T 2 es/Ti es = 1-7 — 2.3. 


7. Discussion 

The HST/COS observations of KOI-81 have revealed the spectrum of the hot companion, 
the stripped down remains of the originally more massive star in the binary. Our derived 
mass, radius, and temperature values for KOI-81 are consistent with previous results. The 
mass estimate of the hot, compact companion is slightly lower than earlier estimates (~ 
0.3Mq; van Kerkwijk et al. 2010), solidifying its place among low mass He WDs and their 
immediate progenitors (with masses below O.SMq where WDs have helium cores that are 
too small to sustain He burning; Silvotti et al. 2012). The radius derived for the subdwarf 
remains essentially unchanged from earlier estimates, and our analysis confirms that it is 
larger than typical for He WDs (see Fig. 15 of Panei et al. 2007) like the other thermally 
bloated WDs found by Kepler (Rappaport et al. 2015). The effective temperature range for 
the hot companion determined through spectral reconstruction is higher than the previous 
estimate of Rowe et al. (2010), placing it in the hot subdwarf (sdB) and WD regime. Of the 
known main sequence (MS) -|- sdB/He WD binaries, KOI-81 has the hottest, low mass WD 
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(the others range between 9-15 kK), althongh several of the sdB/He WD stars in double 
degenerate systems have similar effective temperatures (Brown et ah 2013; Hermes et ah 
2014). 

The relatively long orbital period of KOI-81 suggests that the system widened as the 
result of conservative mass transfer following the mass ratio reversal. The binary probably 
began as a relatively close binary, so that mass transfer was initiated as the donor grew 
in size after concluding core H burning (designated channel 1 in the evolutionary schemes 
presented by Willems & Kolb 2004). Population synthesis models by Willems & Kolb (2004) 
show that the post-RLOF systems may frequently have a remnant mass and orbital period 
like those we hnd for KOI-81 (see their Fig. 2). The progenitor binaries of such systems 
typically have similar and low component masses and initial orbital periods somewhat larger 
than 1 d (see Fig. 2 in Willems & Kolb 2004). For example, van Kerkwijk et al. (2010) 
present an evolutionary scenario for KOI-81 that begins with a 1.8 and 1.3Mq pair with an 
orbital period of 1.3 d. It is curious that the dominant pulsation period is similar to this. 
It may simply be a coincidence, but an alternative possibility is that this pulsation period 
represents the rotational period of the mass gainer in synchronous rotation before the onset 
of RLOF. Mass transfer from the donor would subsequently spin up the gainer, but it might 
take a long time for the angular momentum to be effectively redistributed into the core of 
the gainer star. We speculate that a slower rotating core may be the source of pulsational 
wave generation in the envelope of the gainer B-star. 

The mass donor remnant may have retained enough of it H-envelope to maintain shell 
H-burning longer, so that it evolves to hotter temperature before cooling begins (Maxted 
et al. 2014; Rappaport et al. 2015). Based on the evolutionary scenarios presented by 
Silvotti et al. (2012) and Rappaport et al. (2015), our derived mass of 0.194 ± 0.020Mq for 
the hot companion of KOI-81 is consistent with a low mass, pre-He WD that has not yet 
reached the cooling branch or a He core WD that is already on the cooling branch after 
experiencing episodes of shell H-burning in CNO flashes. However, when plotted in the 
(logTeff, loggf) plane with the evolutionary sequences for low mass white dwarfs of Althaus 
et al. (2013), KOI-81 appears to be in an early cooling phase of He core WD evolution. A 
specihc age estimate is difficult to pin down because the star is located in a mass regime which 
experiences multiple CNO flashes that cause multiple loops in the (log Teg, log gf) diagram. 
Althaus et al. present a grid of masses and cooling ages weighted by the amount of time 
spent in specihc regions of the cooling tracks and the range of masses found in each region. 
Using the values of logTgfj and logg derived by van Kerkwijk et al. (2010), they hnd a mass 
of 0.196 ± O.OOTMq and cooling age of 372 ± 186 Myr for KOI-81 (see their Table 2). Our 
higher values of log Teg ~ 4.3 and logg = 5.8 correspond to a mean evolutionary track mass 
of 0.202 ± 0.015Mq and a cooling age 356 ± 187 Myr, so our derived stellar parameters are 
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broadly consistent with the evolutionary model predictions. 

The evolutionary path of KOI-81 is one of many that create binary systems with low 
mass WDs. The MS + sdB binaries discovered by Kepler (Rappaport et ah 2015) and WASP 
(Maxted et ah 2014) represent those systems with stripped-down remnants formed during 
the first RLOF phase. These occur during the faster part of the pre-He WD evolutionary 
tracks, but it is easier to detect the low mass WDs at this stage because of their relatively 
large luminosity (Istrate et ah 2014). In fact, most of the known low mass He WDs are those 
fainter WDs on the long-lived part of the cooling tracks that are members of binaries with 
even fainter companions. These are often the remanants of a second mass exchange that 
involves a common envelope stage. They are generally short period systems that are sdB -|- 
WD binaries (Heber et ah 2003; Silvotti et al. 2012), but also include sdB -|- sdB systems 
(Ahmad, Jeffery & Fullerton 2004; Kupfer et al. 2015), sdB -|- neutron star companion (Geier 
et al. 2009; Istrate et al. 2014), and sdB -|- substellar companion (Geier et al. 2009; Geier 
2015). 

The detection of the hot companion in KOI-81 provides us with the means to determine 
the stellar parameters in one example of what must be a large population of undetected 
binaries with faint, hot companions (Willems & Kolb 2004; Di Stefano 2011). The fact that 
the B-star component is also a pulsator opens up the possibility to study the stellar interior 
of a star that has been radically changed through mass transfer. Furthermore, an analysis 
of the transit shapes at different phases in the pulsational and rotational cycles may help 
elucidate the nature of the pulsation modes and source of the rotational modulation (Biro & 
Nuspl 2011). The serendipitous discovery of KOI-81 by Kepler has given us the opportunity 
to explore the properties of stars that have survived transformative mass exchange. 
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WAVELENGTH (Angstroms) 

Fig. 1.— Difference spectra of KOI-81 formed by subtracting the occultation phase spectrum 
of the primary alone. The sharp absorption lines from the subdwarf star appear near max¬ 
imum redshift (blueshift) at orbital phase 0 = 0.25 (0 = 0.76). The vertical line segments 
indicate the positions where interstellar lines were removed from the spectrum. 
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Fig. 2.— Top panel; The radial velocity data from Table 1 for the primary star plotted as 
a function of orbital phase. Phase 0.0 corresponds to the time of central transit of the hot 
subdwarf across the face of the brighter B-type star. The asterisk, diamond, and square 
symbols represent measurements from COS, TRES, and KPNO, respectively, and the solid, 
dashed, and dotted lines show preliminary circular hts for the same three sets (with Ki = 7.0, 
6.3, and 8.0 km s“^ and y* = —8.3, —3.4, and —0.5 km s“^, respectively). Lower panel: The 
radial velocity differences (formed by subtracting the systemic velocities for each set reported 
in Table 2) as a function of phase. The shaded region indicates the ±2(T range in the velocity 
curve for the adopted ht of the combined measurements. 
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Fig. 3.— Radial velocity curves for KOI-81 and its companion. The solid circles and open cir¬ 
cles represent the radial velocities derived from COS spectra for the B-star and hot subdwarf, 
respectively. Plus signs represent the measurements of the B-star velocity from ground-based 
spectroscopy (see Fig. 2). 
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Fig. 4.— The reconstructed UV spectrum of the B-star (solid line) depicted in three suc¬ 
cessive wavelength panels. The dashed line shows the COS spectrum obtained when the 
subdwarf was occulted, and it is identical within uncertainties with the reconstructed spec¬ 
trum. A model spectrum constructed from the UVBLUE grid appears as a dotted line above 
the reconstructed spectrum (offset by -1-0.2 in units of normalized flux). Several strong ab¬ 
sorption lines are identihed and marked by the upper vertical line segments, while the lower 
line segments indicate the locations where interstellar lines were removed from the spectrum. 
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Fig. 5.— The reconstructed UV spectrum of the subdwarf star depicted in three successive 
wavelength panels (in the same format as Fig. 4). The top two panels also show for compar¬ 
ison the spectrum of the hot subdwarf CPD—64°481 (solid line offset by -I-1.3 in normalized 
flux). 
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Fig. 6.— The spectral energy distribution of the KOI-81 binary including (from short to 
long wavelength) fluxes from COS (UV), GALEX (NUV), the KPNO 4 m (optical, low 
resolution spectroscopy), 2MASS (near-IR), and WISE (mid-IR). The solid line depicts a 
Kurucz atmosphere flux model for the combined system that is based upon a £t of the 
reddening and the B-star angular diameter. 
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Fig. 7.— Left: The Ca II A3933 feature (structured solid line) from the average of the 
TRES high resolution spectra. A sharp, interstellar component (with a radial velocity of 
Vr = —19.2 ± 0.9 km s“^) is situated near the center of the broad photospheric component. 
The solid and dotted lines show UVBLUE model prohles for Teg = 10.6 and 11.7 kK, 
respectively. Right: The Hy line prohle (solid line) from an average of the KPNO 4 m 
moderate resolution spectra. The other lines represent UVBLUE model spectra for three 
choices of surface gravity. 
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Fig. 8.— The Fourier amplitude spectrum of Kepler short cadence photometry after 
prewhitening the dominant peak at 0.72297 d“^ indicated by the dotted line. The empirical 
noise level is indicated by the red solid line. 
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Frequency (cycles d'^) 


Fig. 9.— Top: The broad differential rotation feature at a frequency of about 1.96 ~ 2.06 
d“^ and the adjacent sharp peak /s = 2.08287 d“^. Bottom: the same diagram but with the 
sharp /s peak prewhitened. 




Fig. 10.— A model representation of the monochromatic intensity in the Aep/er band-pass of 
the rotationally distorted B-star and the small, hot companion star (shown at hrst contact). 
The horizontal gray line shows the derived transit path. 
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Fig. 11.— The phase-folded light curve of KOI-81 for the occultation of the companion 
(plus signs, shifted by half an orbit and offset by -1-0.006 for clarity) and the transit of the 
companion (dots). The model transit light curve is shown by the solid line. 
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Table 1. Radial Velocity Measurements 


Primary / 
Secondary 
(1) 

Date 

(HJD-2,400,000) 

(2) 

Orbital 

Phase 

(3) 

Vr 

(km s“^) 

(4) 

Vr - 7i 
(km s“^) 

(5) 

CT 

(km s“^) 
(6) 

(O-C) 
(km s~^) 

(7) 

Observation 

Source 

(8) 

p 

55136.6338 

0.7248 

0.77 

4.31 

1.57 

-2.34 

TRES 

p 

55139.6337 

0.8504 

2.03 

5.57 

0.85 

0.13 

TRES 

p 

55143.5776 

0.0156 

-5.33 

-1.79 

0.84 

-1.13 

TRES 

p 

55344.8716 

0.4464 

-3.46 

0.08 

3.70 

2.30 

TRES 

p 

55347.8033 

0.5692 

-6.76 

-3.22 

1.43 

-6.05 

TRES 

p 

55350.7684 

0.6934 

3.63 

7.17 

1.62 

0.86 

TRES 

p 

55366.7730 

0.3637 

-8.98 

-7.52 

1.60 

-2.43 

KPNO 

p 

55366.8595 

0.3673 

-11.00 

-7.46 

3.13 

-2.47 

TRES 

p 

55367.8714 

0.4097 

-3.93 

-2.47 

1.59 

1.15 

KPNO 

p 

55369.7871 

0.4899 

-3.17 

0.38 

3.69 

0.80 

TRES 

p 

55374.7716 

0.6987 

2.49 

6.03 

3.23 

-0.36 

TRES 

p 

55375.8748 

0.7449 

6.08 

9.62 

3.44 

2.88 

TRES 

p 

55376.7460 

0.7814 

10.38 

13.92 

2.65 

7.31 

TRES 

p 

55377.8023 

0.8256 

3.48 

7.02 

6.23 

1.03 

TRES 

p 

55380.8825 

0.9546 

-0.89 

2.66 

2.53 

0.76 

TRES 

p 

55401.7202 

0.8274 

5.38 

8.92 

1.86 

2.97 

TRES 

p 

55458.7261 

0.2150 

-7.64 

-4.10 

3.15 

2.47 

TRES 

p 

55459.6061 

0.2518 

-0.91 

2.64 

2.08 

9.37 

TRES 

p 

55469.6640 

0.6731 

4.14 

7.68 

3.54 

1.72 

TRES 

p 

55483.5960 

0.2566 

-16.32 

-12.78 

4.97 

-6.05 

TRES 

p 

55495.5829 

0.7586 

13.85 

17.39 

4.43 

10.67 

TRES 

p 

55734.3079 

0.7571 

-0.33 

8.00 

1.09 

1.27 

HST/COS 

p 

55734.4731 

0.7640 

-1.81 

6.52 

1.09 

-0.19 

HST/COS 

p 

55775.9626 

0.5017 

-9.31 

-0.98 

1.09 

-1.05 

HST/COS 

p 

55841.5269 

0.2478 

-14.29 

-5.96 

1.09 

0.78 

HST/COS 

p 

55841.6811 

0.2542 

-15.88 

-7.55 

1.09 

-0.81 

HST/COS 

p 

56077.7512 

0.1415 

-10.35 

-8.88 

1.84 

-3.65 

KPNO 

p 

56081.7577 

0.3093 

-2.55 

-1.08 

1.87 

5.19 

KPNO 

p 

56082.7568 

0.3512 

-3.22 

-1.76 

2.45 

3.66 

KPNO 

p 

56486.7101 

0.2699 

-17.15 

-15.68 

4.03 

-9.00 

KPNO 

s 

55734.3079 

0.7571 

-102.81 


2.08 

-1.54 

HST/COS 

s 

55734.4731 

0.7640 

-99.37 


2.12 

1.61 

HST/COS 

s 

55841.5269 

0.2478 

101.12 


2.42 

0.15 

HST/COS 

s 

55841.6811 

0.2542 

100.74 


2.81 

-0.21 

HST/COS 
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Table 2. Orbital Elements for KOI-81 


Element 

Value 

P (days) . 

23.8760923“" 

Tt (HJD-2,400,000) 

54976.07186“" 

Ki (km s“^) . 

6.74 ±0.67 

K 2 (km s”’^) . 

101.18 ±0.73 

7 i[TRES] (km s-^) 

-3.54 ±0.62 

7 i[KPNO] (km s"^) 

-1.47 ±0.68 

7 i[COS] (km s-’) .. 

-8.33 ±0.17 

72[C0S] (kms-i) .. 

-0.19 ±0.72 

M 2 /Ml . 

0.0666 ± 0.0066 

asini {Rq) . 

50.91 ±0.47 

rmsi (km s“’) . 

2.9 

rms 2 (km s“^) . 

1.6 


®'Fixed. 


Table 3. Stellar Parameters for KOI-81 


Parameter Primary Secondary 


M/Mq . 2.916 ±0.057 0.194 ±0.020 

R/Rq . 2.447 ± 0.022^^ 0.0911 ± 0.0025'" 

\ogg (cgs) . 4.13*" 5.81*" 

Teff (kK) . 11.7 ±1.5 >19.4 ±2.5 

V sini (km s“^) 296 ±5 <10 


“"Assuming a spherical shape for the primary. 
’"Calculated from M/Mq and R/Rq. 






















Table 4. 


Significant Photometric Frequencies 



Frequency (d 

Amplitude (10 

Phase (rad/27r) 

S/N 

Comment 

/l 

0.722974 ± 0.000002 

257.6 ±4.1 

0.822 ±0.001 

107.4 


/2 

1.32403 ± 0.00001 

33.4 ±2.8 

0.34 ±0.006 

20.6 

k ± k 

fs 

1.08445 ± 0.00002 

22.1 ±3.2 

0.98 ±0.01 

20.5 

fi ± k 

h 

0.96250 ± 0.00001 

34.6 ±3.6 

0.22 ±0.008 

16.7 


fb 

2.08287 ± 0.00001 

31.4 ±2.6 

0.77 ±0.006 

11.4 

frot 

k 

0.36148 ± 0.00002 

47.5 ± 7.7 

0.28 ±0.01 

10.6 

0.5/i 

h 

2.34753 ± 0.00002 

15.7 ±2.3 

0.36 ±0.01 

8.0 


h 

0.70933 ± 0.00009 

6.3 ±4.2 

0.96 ±0.04 

6.9 


k 

1.74647 ±0.00006 

5.8 ±2.3 

0.15 ±0.03 

6.8 


/lO 

2.08357 ± 0.00003 

12.3 ±2.6 

0.37 ±0.05 

4.5 

frot 

/ii 

4.16643 ±0.00008 

4.4 ±2.4 

0.08 ±0.04 

4.3 

‘^frot 

/l2 

0.0837 ± 0.0002 

10.9 ±4.0 

0.25 ±0.02 

3.2 

“^forb 

/l3 

2.08216 ± 0.00003 

10.4 ±2.6 

0.09 ±0.04 

3.1 

frot 






